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SIMPLIFIED THEORY FOR DYNAMIC RELATION OF 
RAMJET PRESSURES AND FUEL FLOW 
By Herbert G. Hurrell 


SUMMARY 

An analysis is made to determine, for control design purposes, the 
approximate response of pressures in the ramjet engine to changes in fuel 
flow. The problem is simplified by separating it into two parts. In the 
first part, the response after dead time is treated by the linearized 
lumped-parameter method. Cons ideration is limited to ramjet operation 
that is critical or supercritical. In the second part, the dead time 
occurring between the change in fuel flow and the beginning of the pres- 
sure response is discussed. 

Comparisons of the theory with experimental data indicate that the 
response after dead time can be approximated by a second-order lead- lag 
and that the dead time can be calculated on the basis of a sonic dis- 
turbance propagating through steady flow. The time constants of the 
theoretical lead- lag depend on the dynamic effects of the diffuser ter- 
minal shock and on the ratios of the significant volumes to the steady- 
state volume flow rate. 


INTRODUCTION 

Recent investigations have indicated that the most suitable controls 
for the supersonic ramjet engine are those that use a pressure, or group 
of pressures, within the engine as the controlled variable. In most ram- 
jets the pressure will be controlled by means of fuel-flow manipulation. 

A method to predict the dynamic relation between these two variables is 
of great interest, therefore, in the synthesis of ramjet control systems. 

When considered in rigorous form, the problem posed by this interest 
is a formidable one. It is concerned with the propagation, reflection, 
and interaction of pressure waves and with the complex flow patterns 
built up by these waves. The variables of such flows axe, of course, 
functions of both time and position, and the part ial differential equa- 
tions relating them are highly nonlinear. Even when a flow can be 


••• •••• ••• 
• + • # 

• • • • • m 

• • # m 

••• ••• 


• •• 
• ••• 


••• ••• 

CONFIDENTIAL 




• • 


• • 

• • • • 
• • 

• • 


• • • 


• •• 




• • 


• • 

• • • 

•• • • 

• •• • 

• • • 


• • • • • • 


• • • 


• • 


: ••• • 
••• •••• 


•••• ••• 


••• 


•• ••• • 

CONFIDENTIAL 




•NAEA RM E57I13 


treated as one dimensional, the nonlinearity of the governing equations 
confronts any attempt to obtain a general solution. In fact, there seems 
to be no convenient way to work with the rigorous problem apart from 
stepwise numerical calculation of each transient. The most convenient 
approach to these numerical calculations is, perhaps, the method of char- 
acteristics, which is well discussed in the literature (e.g., refs. 1 
and 2) . 

A characteristics solution may not be warranted, however, for the 
purpose of synthesizing ramjet controls. Such a solution is very tedious 
and time consuming, and its potential accuracy will not be realized un- 
less the initial conditions of the gas flow are well defined throughout 
the ramjet. Usually such detailed information is not available. Further- 
more, numerical calculations must be interpreted if an understanding of 
the general problem is desired, and unfortunately they tend to hide any 
properties common to all transients. 

With certain stipulations to ensure linearity, of course, the funda- 
mental equations of transient gas flow can be readily integrated (ref. 

3). The equations obtained treat the wave phenomena in approximate form. 
For a realistic treatment of ramjet transients, however, such a solution 
would be very complex and, hence, inconvenient to use in control- system 
design. In addition, interpretation of the solution in terms of changes 
in initial and boundary conditions may again be very difficult. For con- 
trols synthesis, therefore, it appears that further simplification of the 
problem merits consideration. 

The simplification used in the present report is to preclude a wave 
solution altogether. This is done by separating the problem into two 
parts. In the first part, an approximation is sought for the pressure 
response that occurs once the initial disturbance wave has reached the 
sensing station. The approximation is obtained by analyzing the dynamics 
of the ramjet in lumped-parameter fashion, a lumped-parameter analysis 
being one in which the flow variables are treated as functions of time 
alone, rather than as functions of time and position. Such a procedure, 
of course, permits integration of the continuity equation (with respect 
to position) without recourse to Euler’s equation of momentum. In the 
second part, the delay (dead time) between the fuel-flow disturbance and 
the arrival of the initial pressure wave at the sensing station is ana- 
lyzed. This analysis of dead time is made on the basis of sonic propa- 
gation of the wave through steady flow. 

Previous workers in the field of ramjet dynamics, of course, have 
used the lumped-parameter method of analysis (ref. 4, for instance). In 
these works, however, no consideration is given to the movement of the 
diffuser terminal shock that necessarily accompanies any transient con- 
dition of the gas flow. For the present analysis, it was presumed that 
the shock transient makes an important contribution to the dynamics of 
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the ramjet. Accordingly, the position of the shock is treated herein as 
a time -dependent variable. The required definition of the shock tran- * 
sient is obtained from reference 5. 

In the formulation of the simplified theory for response after dead 
time, only critical or supercritical operation of the ramjet during 
transients is treated. The treatment is also restricted to smal i tran- 
sients in order to permit linearization of the equations involved. In 
addition to the presentation of the analysis, this report includes com- 
parisons of the theory with experimental data obtained at the NACA Lewis 
laboratory. Evaluations are made for both parts of the analysis, re- 
sponse after dead time and dead time. 


SYMBOLS 

A cross-sectional area of ramjet 

a speed of sound 

B,C,D steady-state quantities used in eq. (8) 

K gain of pressure in response to fuel flow 

k gain relating change in shock position to downstream pressure 

disturbance 

ky gain relating change in cold volume to change in lumped pressure 

M Mach number 

p static pressure 

R gas constant 

s Laplacian operator 

T static temperature 

t time 

u gas veloc ity 

V volume • 

w gas flow 

Wj, fuel flow 
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x position coordinate in axial direction (+ downstream) 

a dead time 

y ratio of specific heats 

p static density 

a shock time constant 

a' dimensionless shock time constant 

T l> x 2 

time constants 

V T 4 

Subscripts: 

a fuel- injection station 

b location of combustion disturbance 

c cold zone 

e exhaust -nozzle throat 

h hot zone 

ss steady state 

t total conditions 

0 entering shock 

1 leaving shock 


RESPONSE AFTER DEAD TIME 
Lumped-Parameter Analysis 

The lumped-parameter concept of the ramjet can be explained by ref- 
erence to figure 1. As indicated by this sketch, a "hot zone" and a 
"cold zone" are considered, the division being the main flameholding 
element of the combustor. The heat input to the engine is assumed to 
occur only in the hot zone. This zone, which extends downstream to the 
throat of the exhaust nozzle, contains a fixed volume of gas. In the 
cold zone, however, the volume of gas is time -dependent, since the 
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terminal shock of the diffuser forms the upstream boundary. Within each 
zone the variation of gas temperature with position is neglected, a 
value of temperature being considered in each zone. Pressure is lumped 
into a mean value for the combined hot and cold zones. 

In accordance with the lumped treatment of the flow variables, the 
principle of mass conservation states 


w„ = 


a(p c v c : 

dt 


+ v, 


dPfc 
h dt 


(1) 


and if consideration is limited to perturbations from the steady-state 
condition that are small enough to neglect products of perturbation terms, 
the statement reduces to the linear equation 


-Aw_ 


= Pc 


dAV c 

dt 


+ V„ 


<34> c „ ^h 

dt + V h dt 


( 2 ) 


where the A quantities are the small perturbations, and the coefficients 
of the derivatives are steady-state values. (Hereafter, only A quanti- 
ties are time-dependent.) The emission of the term Aw^ from the equa- 
tion results from restricting the analysis to ramjet operation that is 
critical or supercritical during transients. By use of the Laplace trans- 
formation, this linearized continuity statement takes the form 

-Aw e (s) = P c sAV c (s) + V c sAp c (s) + ^sAp^s) (3) 

From the equation of state p = pRT and the assumption that in the 
cold zone the lumped variables of state are related isentropically, the 
changes in density can be related to pressure and temperature changes as 
follows: 


^h( s ) - RT - *2 ^h( s ) 

Ki h RT£ 

(4) 

^P c ( s ) = t rt c 

(5) 


The condition of a choked exhaust nozzle permits the change in gas 
efflux to be related to the pressure and temperature changes in the fol- 
lowing manner: 


Aw e (s) = | Ap(s) - ^ AP h (s) (6) 

In order to define the change in the cold volume, an analysis of 
the shock motion is necessary-. Such an analysis is presented in 
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reference 5, where small perturbations of a normal shock in quasi-one- 
dimensional flow are considered. It is shown in this reference that when 
a disturbance in downstream pressure occurs the resultant change in shock 
position is related to the disturbance by a first-order lag. The time 
constant a and the gain k,™ of this lag are expressed in terms of a 

dimensionless time constant o' as follows: 



, _ r + i a o i 

sp ’ ** a t 0 


(Negative values of k^ denote upstream movement of the shock.) The 

dimensionless time constant a r depends only on the steady-state value 
of the shock Mach number Mq, “being defined in the following manner: 





^2 


1 + 




With this Information from reference 5, the 
ume can be related to the lumped pressure by the 


change in the cold vol- 
equation 


AV C< S > - “v TTVs ^ s) 


(7) 


where 


kv - 

Equations (4) to (7) can now be combined with the continuity expres- 
sion of equation ( 3 ) to obtain a single equation containing only the 
changes in pressure and hot-zone temperature as variables. This equa- 
tion, together with the assumption that hot-zone temperature and fuel 
flow are related quasi- statically, will yield the following transfer 
function for pressure to fuel flow: 

Ap(s) _ K 1 + (2B+o)s + 2Bas 2 ^ 

A w f (6) ! + (B+C+I>fa)s + (B+C )os^ 


••• •• •• 




• * • 

• • • • i 



NACA RM ESttlS 


CONFIDENTIAL 


7 


where 


B 


C 


PV h 


P V c 

r>®T c 


D = 


kyP 2 

wRT c 



A more convenient form of the transfer function is obtained when equa- 
tion (8) is factored as follows: 


Ap(s l C 1 + T l s )( 1 + t 2 s ) 

Aw f (s) (l + t 3 s)(1 + t 4 s) 


where 


T 1 = 0 
t 2 = 2B 


t 3 = 


T, = 


2g(B+C) 


(B+C+D+a) - ^B+C+D+a) 2 * - 4a(B+C) 
2a(B+C) 


(B+C+D+a ) + /y(B+C+D+a) — 4:a(B+C) 

According to lumped-parameter considerations, therefore, the re- 
sponse of pressure to fuel flow in the ramjet is a second-order lead.- lag 


characterized by the time constants t^, t 2 , t^, and t^. These time 

constants, in turn, are functions of the quantities B, C, and D and 
the shock time constant a. The q uan tities B and C can be thought 
of as the "filling times" of the hot and cold volumes, respectively, as 
they define the times required to fill these volumes with gas at the 
steady-state volume flow rate. The quantity D and the time constant 
a represent the contribution of the. shock to the dynamics of the engine. 


These filling times and shock-associated terms can be calculated for 
a given ramjet configuration once the steady-state flight conditions and 
engine operating point are specified. The simplified description of the 
steady flow necessary for these calculations is readily obtainable; 
either the well established theory for steady flow or experimental data 
can be used for this purpose. For computing the shock time constant a 
(as well as the gain ky. involved in the quantity D) the curve shown 
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in figure 2 is useful. This curve, which is reproduced from reference 5, 
gives the value of the dimensionless shock time constant o' for shock 
Mach numbers ranging in value from 1.0 to 3.0 (y* = 7/5). 

It should he emphasized that the transfer function derived in this 
section must he combined with a term for dead time to obtain the complete 
transfer function for the ramjet. The complete transfer function is 

Ap ( s ) = Ke -as (1 + TpsKl + t 2 s) 

Aw f ( s ) (1 + t 3 s)(1 + 

where a is the dead time at the station of interest. The discussion 
of dead time follows the experimental corroboration of the foregoing 
analysis . 


Comparison of Theoretical and Experimental Responses 

The lumped-parameter theory for response after dead time has been 
compared with the experimentally determined responses of two ramjets. 
These engines were both full-scale and representative of current ramjet 
design. They were considerably different, though, in size, configura- 
tion, and flight specifications and thus afforded a rather extensive 
check on the theory. The dynamic behavior of these engines had been 
well defined experimentally. Test programs had been conducted in which 
steps or sinusoidal variations in fuel flow were imposed by means of 
specially designed electrohydraulic servo systems. The resultant re- 
sponse of pressure was measured at several engine locations with instru- 
mentation suitable for recording pressure fluctuations up to 100 cycles 
per second in frequency. During the tests, the ramjets were operated in 
supersonic airstreams generated by free- jet facilities. 

In the calculations of the theoretical responses, the steady-state 
descriptions required were obtained by means of one-dimensional theory 
in conjunction with the available experimental data. For simplicity, 
the assumption was made that the entire heat input to an engine occurred 
at its main flameholding element. In the computations of the lumped val- 
ues of flow variables, the variables were first determined as functions 
of position in the region being considered, and then the means were taken 
by a volume-weighting process . In order to determine the values of the 
shock time constant a and the gain k^ for one engine, designated 

ramjet I, the parameters involved in these quantities were averaged over 
the region traversed by the shock. This was necessary because rather 
large movements of the shock were being treated. For the other engine, 
designated ramjet II, the shock terms could be computed from the initial 
condition, since the shock travel was small. 
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The responses calculated for the two engines are compared with ex- 
perimental data in figure 3. The comparison, which is presented in fig- 
ures 3(a) and (b) for ramjet I and in figure 3(c) for ramjet II, is made 
on the basis of frequency response; amplitude ratio and phase shift are 
shown as functions of frequency. The amplitude ratios are given in nor- 
malized form, that is, they have been adjusted to a steady-state gain of 
unity. The responses from the lumped-parameter theory are shown in the 
figure by the solid curves, while the experimental data are represented 
by symbols. The straight-line approximations of the amplitude curves 
are given as dashed lines, and the break frequencies (l/2«T ) associated 
with the calculated time constants are indicated. 

From these break frequencies it is evident that the dynamics of the 
shock play an important ' role in the theory formulated. For instance, 
the time constant r^, which is equivalent to the shock time constant a, 

contributes strongly to the shape of the theoretical amplitude curve for 
each ramjet . 

The experimentally determined responses include the effects of dead 
time, of course. Since the amplitude component of dead time is always 
unity, the amplitude data afford a direct means of evaluating the theory 
for response after dead time. In plotting the phase-shift curves, how- 
ever, it was necessary to combine the phase characteristics of the lumped- 
parameter theory with those of the appropriate dead times. The values of 
dead time used for this purpose were determ in ed in the step-response 
tests . The phase curves without this adjustment are also included in the 
figures . 

The data for ramjet I in figure 3(a) are presented for one location 
of the pressure tap but for fuel-flow sinusoids of two different ampli- 
tudes . The data shown by the circular symbols were obtained with the 
fuel-flow amplitude as small as practical; the square symbols represent 
data from a sinusoid of large amplitude. The figure shows the experi- 
mental results for both sinusoids are described reasonably well by the 
theoretical curves. In fact, the change in fuel-flow amplitude produced 
no sizable effect on either the amplitude or phase characteristics of 
the pressure response. On the basis of this figure, therefore, the lin- 
earization performed in the analysis by the assumption of small transients 
does not impair the usefulness of the theory. The theory may be used to 
predict the approximate response for large as well as small transients. 

Any difference in response from one station to another is not con- 
sidered, of course, in the lumped-parameter analysis. The extent of 
error in such a treatment can be seen in figure 3(b). The shaded band 
in the figure indicates the effect of pressure-sensor location an the 
experimental results for ramjet I. This effect is shown only for the 
amplitude characteristics because the phase effect was predominantly due 
to dead-time changes. The width of the shaded band at a given frequency 
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denotes the maximum variation of the amplitude -ratio data for stations 
throughout the ramjet. (The variation appeared to be completely random.) 
The figure indicates that the theory provides a reasonable approximation 
to the response at any station in the ramjet. Further evidence of this 
is shown in figure 3(c), where the amplitude data for ramjet II are pre- 
sented for two stations in the engine. In this figure the curve for the 
theoretical response gives a good approximation of the experimental data 
for both stations over most of the frequency range. 

In summary of the experimental corroboration, figure 3 indicates 
that the lumped-parameter theory reasonably approximates the response of 
ramjet pressures to fuel flow, exclusive of dead time, of course. Evi- 
dence is given that the theory may be used for pressures throughout the 
engine and for large as well as small disturbances in fuel flow. It is 
believed that the theory is sufficiently accurate to be useful in the 
design of ramjet controls. 


DEAD TIME 

The dead time in the response of ramjet pressures to fuel flow con- 
sists of two delays: first, the time required for the fuel (injected at 

a new rate) to travel from the injection station to its place of burning 
in the combustion chamber, and second, the time spent by the resultant 
pressure wave in propagating from the combustion disturbance to the 
pressure-sensing station. There appears to be no way to locate, by ana- 
lytical means, the place in the combustion chamber where combustion is 
disturbed and the pressure wave originates. Experimental data indicate 
the location to be some distance downstream of the main flameholder; 
however, the distance depends on the initial state of the gas flow, the 
size of the fuel-flow change, and whether the fuel flow was increased or 
decreased. Fortunately, the error caused by this uncertainty is usually 
small in comparison to the dead times being calculated. Therefore, an 
approximate location may be used. For ramjets I and II the distance was 
assumed to be one-third of the way from the main flameholder to the noz- 
zle throat. 

Once an approximate location has been selected for the combustion 
disturbance, the calculation of dead time becomes essentially a problem 
of determining the gas velocities and the sonic velocities in the engine 
for steady- state operation. The fuel is considered to travel at the ve- 
locities of the gas, and the pressure wave is assumed to propagate at 
sonic velocities relative to the gas. The assumption that these veloc- 
ities are those of the steady state is, in general, a small-perturbation 
treatment; for the general problem relating to control response, of 
course, the fuel-flow change may be preceeded by transient flow 
conditions . 
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Dead time may be defined in terms of the steady- state values of the 
gas velocity u and sonic velocity a as follows: 



where a x is the dead time at station x, Xg is the fuel injection 
station, and x^ is the station of combustion disturbance. The veloc- 
ity u is the component in the x-direction, and the value of x is con- 
sidered to increase in the downstream direction. The plus signs apply 
when x is greater than x^, and the minus signs are to be used when x 
is less than x^. 

Methods for calculating the velocities u and a for steady flow 
are well known and need not be discussed herein. The tables of refer- 
ence 6 are very convenient for this purpose. Often the theory for these 
calculations can be supplemented by experimental data. This was done 
in the velocity calculations for ramjets I and II; seme experimental 
pressure and temperature data were used in conjunction with one- 
dimensional theory. In addition, an in the calculations for the lumped- 
parameter response, it was assigned for simplicity that steady-state com- 
bustion occurred only at the main f lameholders . 


The dead times calculated for ramjets I and II are compared in fig- 
ure 4 with those determined experimentally during the step-response 
tests. The dead times are given for the region of each engine downstream 
of the normal shock. The calculated variation of dead time in this re- 
gion is shown by the curve, and the experimental data are represented by 
the symbols. For ramjet II the data are given for both a step increase 
and a step decrease in fuel flow. For ramjet I there was no separation 
of the data with the direct ion of the step . 

As can be seen in figure 4, the calculated dead times agree reason- 
ably well, in general, with the experimental values. Some deviation, 
however, can be seen for ramjet II (fig. 4(b)). For this engine, in 
which the data depend on the direction of the fuel-flow step, the cal- 
culated curve agrees better with the dead times measured for a step in- 
crease than with those measured for a step decrease. The deviation of 
the curve from the step-decrease data is most pronounced for stations 
downstream of the flameholder. For stations upstream of the flameholder, 
the curve and the step-decrease data are in fair agreement . 
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CONCLUDING REMARKS 

For critical or supercritical operation of the ramjet, it appears 
that disturbances in fuel flow provoke pressure responses that, exclusive 
of dead time, can be described approximately by a second-order lead- lag 
function obtained from lumped-parameter considerations. The approxima- 
tion is considered sufficiently accurate to be useful to the controls 
designer. The time constants of the theoretical lead- lag are completely 
determined by the geometry of the ramjet and a simplified description of 
the steady-flow conditions preceeding the transient. These time con- 
stants, therefore, can be readily calculated with the use of steady-state 
data or well established flow theory. 

The time constants are composed of terms that describe the filling 
times of the ramjet and the dynamic effects of the diffuser terminal 
shock. Filling time refers to the time required to fill a significant 

volume in the engine at the volume flow rate that exists in the steady 

state. The dynamic effects of the shock result from the transient in 
shock position that accompanies the pressure response, the shock position 
being dynamically related to the pressure by a first-order lag. 

The complete transfer function relating a ramjet pressure p to 
fuel flow wp is obtained by multiplying the second-order lead- lag 
function by the term for the dead time at the station involved. The 
complete transfer function is 

Ap(s) _ „ -a U 1 + T l s ) ( 1 + t 2 s ) 

Aw f (s) e (l + T 3 s) (l + T 4 s) 

where a is the dead time, \, ^ 3 , and are time constants, K 

is the gain, and s is the Laplacian operator. The dead time can be 
calculated once the place of the combustion disturbance has been located 
with reasonable accuracy. The fuel can be considered to travel from the 
injector to this place of disturbance at the steady-state gas velocities; 
the resultant pressure wave can be assumed to propagate relative to the 
gas at the steady-state sonic velocities. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 16, 1957 
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Figure 1. - Lumped- parameter concept of ramjet. 
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# • Figure 2. - Variation of dimensionless shock time constant with shock Mach number. 

••••• Ratio of specific heats, 1.4. 
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Normalized amplitude ratio 
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(b) Ramjet I; range of pressure tap locations. 

Figure 3. - Continued. Comparison of theoretical and experimental responses of pressure to fuel flow after dead 
time. Diffuser total- pressure recovery , 95 percent of critical. 
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Frequency, cps 


(c) Ramjet II. Amplitude of fuel-flow sinusoid, 12.5 percent of critical fuel flow. 

Figure 3. - Concluded. Comparison of theoretical and experimental responses of pressure to fuel flow after dead time. 
Diffuser total-pressure recovery, 95 percent of critical. 
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(a) Ramjet I. 
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(b) Ramjet II. 

4. - Comparison of calculated and experimental dead times of pressure to fuel 
. Diffuser total-pressure recovery, 98 percent of critical; step change in 
flow, 10 percent of critical fuel flow. 
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